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Abstract 

Dynamic relaxation of the XY model quenched from a high tempera- 
ture state to the critical temperature or below is investigated with Monte 
Carlo methods. When a non-zero initial magnetization is given, in the 
short-time regime of the dynamic evolution the critical initial increase of 
the magnetization is observed. The dynamic exponent is directly deter- 
mined. The results show that the exponent varies with respect to the 
temperature. Furthermore, it is demonstrated that this initial increase of 
the magnetization is universal, i.e. independent of the microscopic details 
of the initial configurations and the algorithms. 
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1 Introduction 



Recently much attention has been drawn to the short-time universal scaling be- 
haviour of critical dynamic systems. A typical example is the dynamic relaxation 
of magnetic systems quenched from a high temperature state to the critical tem- 
perature. After a time period enough long in the microscopic sense, during which 
the non-universal short wave behaviour is swept away, the universal scaling be- 
haviour appears. Such a time period is called the microscopic time scale tmic- 
tmic is in general very small compared with the macroscopic time scale, which is 
typically tmac ~ t"''^ or tmac ~ L^- Here r is the reduced temperature and L 
is the lattice size of the systems, while u and z are the critical exponents. At 
this macroscopic early stage of the time evolution, important is that even though 
the spatial correlation length of the systems is still very small in the macroscopic 
sense, due to the infinite time correlation length the dynamic systems already 
evolve universally. This is out of the traditional belief that the universal dy- 
namic scaling emerges only when the spatial correlation length becomes very 
big. 

One interesting phenomenon is that, if a non-zero magnetization is given to 
the initial state at very high temperature, at the macroscopic early time the 
magnetization surprisingly undergoes a critical initial increase [|l], § 

M(t)~mot^ (1) 

where ^ is a new critical exponent which is independent of the static exponents (3, 
V and the dynamic exponent z. The exponent Q is related to the dimension xq of 
the initial magnetization mg by ^ = (xq — (3/v)/z. The magnetization continues 
the increase in a time scale t ~ itlq , then reaches its maximum and crosses 
over to the well known long-time universal behaviour. 

Numerically the critical increase of the magnetization has directly been ob- 
served 1^, H, 0, ||. The critical exponent 6 was determined to a satisfactory 
accuracy for the two dimensional Ising model and Potts model as well as the 
6-state clock model [0]. Numerical results show rather clean short-time scaling 
behaviour. However, all these models are relatively simple and the spins locate 
in discretized configuration spaces. It would be very interesting to investigate 
whether there exists short-time scaling behaviour in more complex models. 

More important, up to now almost all the numerical simulations for the short- 
time dynamics have been done with the heat-bath algorithm. For the heat- 
bath algorithm the microscopic time scale tmic for the two dimensional Ising and 
Potts model is not bigger than one or two Monte Carlo time steps. This is a 
miracle. If one Monte Carlo time step is typically the microscopic time unit, one 
would expect that tmic should be around 10 ~ 50 Monte Carlo time steps. Even 
though the scaling form in the short-time dynamics and its applications have 
extensively been investigated, it is already overdue to understand universality in 
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the short-time dynamics, i.e. whether the short-time scahng behaviour is really 
independent of the algorithms, lattice types and other microscopic details. 

In a recent paper [^, the authors have carefully analysed the short-time 
behaviour of the critical dynamics for the two dimensional Ising model and 
Potts model. Simulations have been performed with both the heat-bath and 
the Metropolis algorithm. Indeed, it is only by chance that for the heat-bath al- 
gorithm tmic is negligibly small. Actually for the Metropolis algorithm, tmic ~ 30. 
Within the time period up to tmic, the dynamic behaviour for the Metropolis al- 
gorithm is very different from that of the heat-bath algorithm. When t > tmic, 
however, the dynamic systems with both algorithms present the same universal 
behaviour. The measured values of the critical exponent 6 are compatible within 
the statistical errors. This is a first step to the verification of universality in the 
short-time dynamics. More understanding is urgent and important 

In this paper, we will report results for the numerical simulations of the short- 
time dynamics of the two dimensional XY model. Special attention will be put 
on universality. In the next section a short description of the XY model is given. 
In the section 3, the numerical data are presented and conclusions are given in 
section 4. 



2 The XY model 

The XY model in two dimensions is defined by the Hamiltonian 

H = KY.S.-S,, (2) 

<ij> 

where Si = {Si^^, Si^y) is a planar unit vector at site i and the sum is over the 
nearest neighbours. In our notation the inverse temperature has been absorbed 
in the coupling K. 

The XY model is the simplest statistical system which exhibits a continu- 
ous symmetry. It is known that at a certain critical temperature the XY model 
undergoes a Kosterlitz-Thouless phase transition |Tl|, [1^. Near the critical tem- 



perature the spatial correlation length diverges exponentially rather than by a 
power law as in the normal second order phase transition. Below the critical 
temperature the system remains critical in the sense that the spatial correlation 
length is divergent. No real long range order emerges in the XY model. The 0(2) 
symmetry of the XY model stays unbroken in the whole temperature regime. 

The XY model is a very important model since it describes the critical prop- 
erties of the superfluid hehum. It is closely related to the 0(2) cr-model in field 

"'^Very recently some discussions on universality in short-time dynamics for the two dimen- 
sional Ising model have also been made with respect to the lattice types and update schemes 
even though the critical exponent 9 there was not confidently extracted due to relative small 
lattices or some other reasons [|[ p^ . 
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theory. Its generalization such as fully frustrated XY model attracts more and 
more attention |]13l O, Il5[. The XY model is also a good laboratory to study 



the more general Heisenberg model. 

However, due to the exponentially divergent spatial correlation length, the 
numerical simulation of the XY model is very difficult. Even though some papers 
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can be found concerning the dynamic properties of the XY model 
knowledge about the dynamic XY model is still very poor. In this paper we will 
investigate the short-time universal behaviour of the dynamic XY model. We will 
concentrate on the critical initial increase of the magnetization and measure the 
critical exponent 6. Special attention will be put on universality. Different from 
the case of the Ising model, where the spins take only the values ±1, the spin 
configuration space for the XY model is a unit circle. This non-trivial configu- 
ration space allows us to investigate whether the short-time universal behaviour 
depends on the microscopic details of the initial configurations. 



3 Numerical simulations 

Following Janssen, Schaub and Schmittmann's idea [jll , we investigate a dynamic 
relaxation process starting from an initial state with very high temperature and 
small magnetization. As a first approach to the dynamic XY model, we do not 
consider the effects of the vortices. The very high initial temperature requests 
that the spin at each lattice site is generated independently. However, the way 
how to generate a non-zero initial magnetization in a certain direction is not 
unique. A natural way is to introduce an initial external field, e.g. in the x 
direction, as it was used in the numerical simulation of the clock model Then 
the initial Hamiltonian, i.e. that for generating the initial magnetization, can be 
written as 

Hoi=2hY,S,,,. (3) 

i 

If we define the magnetization as 

M{t) = ^,J:S^ (4) 

with L being the lattice size, the initial Hamiltonian Hqi gives an initial magne- 
tization 

M(0) = (mo, 0) ^ (/i, 0), 0. (5) 

In this paper we are only interested in the case of small mo. 

To prepare the initial state, we update the system described by the initial 
Hamiltonian Hqi until it reaches equilibrium. Then the generated configurations 
of this initial system are used as the initial configurations of the dynamic system. 
If the lattice size is infinity, in each initial configuration an exact value (mo, 0) of 
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the initial magnetization M(0) is automatically achieved. However, the practical 
lattice size is finite and the initial magnetization M(0) fiuctuates around (mo, 0). 
This is a kind of extra finite size effect. It causes a problem in a high precision 
measurement. In order to reduce this effect, a sharp preparation technique has 
been introduced to adjust the initial magnetization in the numerical simulations 



of the Ising model and the Potts model ^, |^, [T^ : we randomly take one spin 
on the lattice and fiip the spin if the updated magnetization comes nearer to the 
expected value; we repeat this procedure until the expected initial magnetization 
is achieved. Numerical data show that the sharp preparation technique improves 
efficiently the results and especially helps to obtain better results in relatively 
small lattices. 

In the numerical simulation of the XY model, we also implement the sharp 
preparation technique. However, due to the fact that spins in the XY model 
are planar unit vectors, this procedure becomes slightly more complicated. We 
proceed in the following way: 

(i) If the configuration generated by the initial Hamiltonian Hqi does give 
the value < Sx >= m but not mo, we update a randomly chosen spin. If the 
resulting magnetization < Sx > is nearer to mo, we accept it otherwise keep the 
old configuration. We continue in this way until the difference |m — mo| < S with 
6 being a certain given small value. In our simulations we take 6 to be 2.5 percent 
of mo; 

(ii) After having adjusted < Sx >, we turn to the magnetization < Sy >. If 
I < 5^, > I > 5, we randomly select a spin Si = {Si^x, Si^y) and change the sign of 
Si^y. If I < Si^y > I becomes smaller, we accept the new configuration otherwise 
keep the old one. We continue until \ < Sy > \ < 6. In this procedure (ii), the 
value < Sx >= m already prepared in (i) remains unchanged. 

After the preparation of the initial configuration, the system is released to 
a dynamic evolution at the critical temperature or below with the Metropolis 
or the heat-bath algorithm. We have performed the simulations for lattice sizes 
L = 8, 16, 32, 64 and 128. The magnetization is measured up to Monte Carlo 
time step t = 150. The average is taken over 40 000 samples with independent 
initial configurations for the lattice size L < 64 and 12 000 — 30 000 samples for 
the lattice size L = 128 depending on the initial magnetization. For smaller mo 
we take relatively large statistics. Errors are estimated by dividing the data into 
three or four groups. In this paper we take the critical temperature of the XY 
model from the literature |TP[, Tc = 1/Kc = 0.90. Unless we explicitly specify, 
all the discussions below are assumed to be at the critical temperature. 

In Fig. m, the time evolution of the magnetization at the critical temperature 
with the initial magnetization mo = 0.02 for the Metropolis algorithm is displayed 
for different lattice sizes. In the figure M{t) is the x component of the magne- 
tization M{t). The y component of the magnetization M{t) remains zero since 
the initial value is zero. From the figure we can see that for L = 64 the finite size 
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effect is already very small and the curve almost completely overlaps with that 
of L = 128. We have discussed above that the universal behaviour appears only 
after a microscopic time scale t^ic- Theoretically one would expect and it is also 
observed in the cases of the two dimensional Ising model and Potts model that 
tmic is in general around 10 ~ 50. In Fig. |l] one can see this clearly. In the first 
30 time steps there is no universal power law behaviour but after that it indeed 
appears. From the slope of the curve one may measure the critical exponent 6. 

How important is the sharp preparation of the initial magnetization? This 
practically depends on how big the lattice size and how small the initial magne- 
tization mo is. On the other hand, since the exponent 6 is defined in the limit 
mo = 0, the practically measured exponent 6 from the power law behaviour (|l|) 
shows in general some weak dependence on mo when mo is finite. The stronger 
the dependence of 6 on mo is, the more important becomes the sharp preparation. 
In Fig. the magnetization without a sharply prepared initial magnetization is 
displayed for different lattice sizes. For comparison, the dotted line shows that 
with a sharply prepared initial magnetization for the lattice size L = 64. Com- 
paring Fig. 1^ and Fig ^ we see that the difference between the curves with and 
without the sharp preparation of the initial magnetization becomes already quite 
small when the lattice size reaches L = 64. Such a small difference is also partly 
due to the quite weak dependence of 6 on mo, which can be seen later. 

With the sharp preparation of the initial magnetization, the exponent 6 mea- 
sured from lattice size L = 64 and 128 are 6 = 0.250(1) and 0.249(4) respectively. 
Within the statistical errors we already can not distinguish the results for the 
lattice size L = 64 and L = 128. Without the sharp preparation of the initial 
magnetization, we get the exponent 6 = 0.252(2) for lattice size L = 64, which 
shows a slightly bigger value and fluctuation compared with that with the sharp 
preparation of the initial magnetization even though the difference is small. In 
the following simulations, the sharp preparation technique is always adopted. 

Is the power law scaling behaviour (|I|) really universal? Would it depend 
on the microscopic details of the initial configurations, algorithms and lattice 
types and so on? In this paper we will show that the power law behaviour is 
indeed independent of the microscopic details of the initial configurations and 
the algorithms. 

In order to generate an initial state with a non-zero magnetization, using the 
initial Hamiltonian Hqi given in (^) is a natural way but by no means unique. 
An example of alternative methods may be the following: in each lattice site, 
the spin orients towards the pure positive x direction {Si^^ = 1, Si^y = 0) with 
a certain probability, otherwise randomly. This initial state can be described by 
an initial Hamiltonian 



Here the angle 0j is defined by Si^x = cos (pi and Si^y = sin (pi. Properly choosing 
the constant C2 one obtains the expected initial magnetization mo. Another 




(6) 
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possibility is: assuming that the orientation of initial spins is restricted to either 
pure X or pure y direction, we give a slightly bigger probability to generate spins 
in the positive x direction than in others. The corresponding initial Hamiltonian 
is 

Hoz = E 111 (C35(0^) + <5(0^ - vr) + 5(0, - 7r/2) + 5(0, + 7r/2)). (7) 

i 

It is clear that the preparation of the initial configurations given by Hq2 and i^oa 
is rather simple [P, ||. 





Metropolis 


mo 


0.02 


0.01 












e 


.249(4) 


.252(4) 


.248(4) 


.252(7) 





Heatbath 


mo 


0.01 




Hqi 


e 


.253(5) 



Table 1: The exponent 6 measured for lattice size L = 128 with different types 
of initial configurations and algorithms. 



T 


0.90 


0.86 


0.70 


0.50 


0.30 


9 


0.250(1) 


0.264(5) 


0.287(3) 


0.283(4) 


0.282(1) 



Table 2: The exponent 6 measured for different temperatures with the Metropo- 
lis algorithm. The lattice size is L = 64. 

In Fig. |] the time dependent magnetization for different types of initial con- 
figurations is plotted for the lattice size L = 128. The solid lines above and below 
are the results for an initial magnetization mo = 0.02 and mo = 0.01 generated 
from Hqi. The dotted line is the magnetization with mo = 0.02 from H02, and 
the dashed line corresponds to that of mo = 0.01 with ifos- In Table 1 the corre- 
sponding 6 measured in a time interval [40, 150] are listed. We see that all three 
initial Hamiltonians give almost the same results. The difference of the initial 
configurations is swept away in more or less one Monte Carlo time step. Fur- 
thermore, the difference of 6 measured from different initial magnetizations mo is 
quite small and already within the statistical errors. Therefore the extrapolation 
of 6 to the limit mo = is not necessary here. This is also one of the reasons why 
the results with and without the sharp preparation of the initial magnetization 
are not so different. 

Before we continue the discussions of the numerical data, we would like to 
make some comments here. For years it is believed that two exponents (3 and u 
sufficiently describe the critical scaling properties of most magnetic systems in 
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equilibrium and the dynamic scaling properties can be described by the dynamic 
exponent z. An essential point in the short-time dynamic scaling exists in the 
claim that an independent critical exponent Xq (or 9) should be introduced to 
specify the dependence of the scaling behaviour on the initial magnetization. In 
principle, however, there are other choices for the scaling variable. For example, 
the initial magnetic field h in Hqi may also be used. Some discussions concerning 



what is a better choice of the scaling variable has recently been made P0| , [18 
5TI]. In our numerical simulations, we have demonstrated that a non-zero initial 
magnetization mo can be realized in different ways, by introducing either an initial 
magnetic field with Hqi or some other initial systems described by Hq2 oi Hq^. 
However, the exponent 6 is the same for different types of initial configurations. 
Therefore these different ways may be considered as the microscopic details for 
the initial state. Introducing an initial magnetic field is only one possibility to 
generate a non-zero initial magnetization. In this sense, the scaling variable mo 
seems to be more general. 

Now let us come back to our discussions of the numerical results. In Fig. § the 
time evolution of the magnetization with mo = 0.01 with the heat-bath algorithm 
is also displayed by the solid line in between. In the first 20 ~ 30 time steps, 
its behaviour is different from that with the Metropolis algorithm. After that, 
however, as in the case of the Metropolis algorithm it stabilizes to the universal 
power law behaviour. To see this more clearly, in Fig. § we have plotted the 
exponent ^ as a function of the time t for both the heat-bath and the Metropolis 
algorithm. The exponent 6 at time t is measured as the slope of the curve in a 
time interval [t,t + 20]. Error bars are estimated by dividing the total sample 
into three groups. After a microscopic time scale tmic ~ 20 — 30, the exponents 
6 for both the heat-bath and the Metropolis algorithm overlap each other. The 
relatively small error bars for the exponent 6 at certain time periods may come 
from the fact that the errors are estimated from only three groups of the data. 
The final values for 6 are given in Table 1. The results for both the heat-bath 
and the Metropolis algorithm are also consistent within the statistical errors. All 
these results strongly support universality in the short-time dynamics. 

Finally we have also performed the simulations with the Metropolis algorithm 
for the temperature below the critical temperature. Since the XY model remains 
critical, a similar scaling form is expected. In Fig. |^, the magnetization for L = 64 
and different temperatures is plotted versus time t in double-log scale. For the 
temperatures T = 0.90 and 0.86 the initial magnetization is mo = 0.02. For the 
temperatures T = 0.70, 0.50 and 0.30 the initial magnetization is mo = 0.01. As 
before, the weak dependence of the exponent 6 on mo has not been considered 
since it is within our statistical errors. The exponent 6 measured in a time interval 
[40, 150] for different temperatures is listed in Table 2. It is known that in the 
equilibrium the critical exponents in general depend on the temperature. Here 
we see the exponent 6 also varies with respect to the temperature. This situation 
is similar to that of the clock model 0. 
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4 Conclusions 



We have numerically investigated the short-time behaviour of the dynamic relax- 
ation of the two dimensional XY model at the critical temperature and below, 
starting from an initial state with a very high temperature and non-zero mag- 
netization. The critical initial increase of the magnetization is observed and 
the exponent 6 is determined. The results show that as the temperature de- 
creases, the exponent 9 first increases rather rapidly and then decreases slowly. 
The independence of the scaling behaviour on the microscopic details of the ini- 
tial configurations and the algorithms is demonstrated and the microscopic time 
scale tmic ~ 30. Universahty in the short-time dynamics is confirmed. Further 
extension of this work remains important, such as the determination of the crit- 
ical temperature and the static exponents from the short-time dynamics and an 
investigation of the effects of the vortices. 
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Figure 1: The time evolution of the magnetization for the XY model with 
mo = 0.02 for different lattice sizes with the Metropolis algorithm is plotted in 
double-log scale. M{t) is the x component of the magnetization M{t). The sharp 
preparation technique for the initial magnetization is adopted. 
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Figure 2: The time evolution of the magnetization for the XY model with 
mo = 0.02 for different lattice sizes with the Metropolis algorithm is plotted in 
double-log scale. The sharp preparation technique for the initial magnetization is 
not adopted. The dotted line shows the magnetization with the sharp preparation 
for L = 64 for comparison. 
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Figure 3: The time evolution of the magnetization for the XY model for lattice 
size L = 128 with different types of initial configurations and algorithms is plotted 
in double-log scale. The sharp preparation technique for the initial magnetiza- 
tion is adopted. The solid lines are those obtained with the initial Hamiltonian 
Hoi. The sohd hnes above and below correspond to mo = 0.02 and 0.01 with 
the Metropolis algorithm. The solid line in between is that for the heat-bath 
algorithm with niQ = 0.01. The dotted line is for the Metropolis algorithm with 
mo = 0.02 prepared with Hq2 and the dashed line is for the Metropolis algorithm 
with mo = 0.01 prepared with ifos- 
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Figure 4: The exponent 6 measured as a function of time t for both the heat- 
bath and the Metropohs algorithm with tjiq = 0.01 and L = 128. 6{t) is obtained 
in a time interval [t,t + 20]. The sharp preparation technique for the initial 
magnetization is adopted. 
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Figure 5: The time evolution of the magnetization for lattice size L = 64 and 
different temperatures with the Metropolis algorithm is plotted in double-log 
scale. The sharp preparation technique for the initial magnetization is adopted. 
The temperature parameters are T — 0.90, 0.86, 0.70, 0.50 and 0.30 (from below). 



14 



